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Abstract 

Pairwise forces between particles in cosmological N-body simulations ai^e generally soft- 
ened to avoid hard collisions. Physically, this softening corresponds to treating the particles 
as diffuse clouds rather than point masses. For particles of unequal mass (and hence un- 
equal softening length), computing the softened force involves a nontrivial double integral 
over the volumes of the two particles. We show that Plummer force softening is consistent 
with this interpretation of softening while spline softening is not. We provide closed-form 
expressions and numerical implementation for pairwise gravitational force laws for pairs 
of particles of general softening scales ei and £2 assuming the commonly used cloud pro- 
files: NGP, CIC, TSC, and PQS. Similarly, we generalize Plummer force law into pairs of 
particles of general softenings. We relate our expressions to the gaussian, Plummer and 
spline force softenings known from literature. Our expressions allow possible inclusions of 
pointlike particles such as stars or supermassive black holes. 
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1 Introduction 



Gravitational softening is a building block in the foundation of any modem cos- 
mological N-body simulation; at the same time it has traditionally been given rela- 
tively little coverage in literature. In this paper we give this problem a share of our 
dedicated attention. 



Any practical numerical simulation will have finite dynamic range, limiting the res- 
olution with which the system of interest may be studied. Oftentimes, the resolution 
requirements for simulations are not completely uniform across the simulation vol- 
ume: there may be special regions of interest, which we would like to study at high 
resolution, whereas the remainder of the simulation volume need not be simulated 
in great detail. For example, we might be interested in simulating individual dark 



Preprint submitted to Elsevier 



9 April 2008 



matter halos embedded within realistic, time-varying cosmological environments 

(e.g. m, m). 



Given finite computational resources, an efficient strategy is to employ adaptive 
resolution, i.e. to use high resolution in the regions of interest and low resolution 
elsewhere |l3l. Two types of resolution that arise in N-body simulations are mass 
resolution and force resolution. The mass resolution is simply limited by the finite 
number of particles; increasing the particle number within some region increases its 
mass resolution. Force resolution is limited by softening of pairwise forces between 
particles. For simulations of coUisionless systems, pairwise forces must be softened 
at short distances to avoid artificial colli sionality due to finite particle number (e.g. 
[IH, [El, [|6l and references therein). For example, Plummer softening corresponds 
to calculating the force between two particles separated by r as 

-> mim2 



where ep here is the Plummer softening length. Other softening laws are commonly 
used in the literature, with their own corresponding softening lengths. Increasing 
the force resolution therefore means decreasing the force softening length. This is 
allowed within regions of high mass resolution where the particle number is en- 
hanced; elsewhere, the force softening length must remain large to avoid coUision- 
ality. One common choice is to scale the softening length as the cube root of the 
mass, e oc m^^^, which holds fixed the maximal density of all particles. 

In Eqn. ([U), one subtlety that arises when computing the pairwise force between 
particles of different mass is that the appropriate choice of softening length be- 
comes uncertain: should the force be softened by the larger e, the smaller, or some 
average of the two? The chief requirements are that the pairwise forces be sym- 
metric in order to conserve momentum (i.e. the same e is used for F12 as for F21); 
and that the softening length is not too small, in order to avoid hard collisions. For 
example, taking e = max(£:i, 62) satisfies both of these conditions. This choice is 
clearly not optimal, however, since the effective force resolution is degraded. 

In addition to basic conservation laws, the consistency of computed forces with 
gravity is also important. Consider a set of point masses connected by massless 
springs for example: such system clearly satisfies both energy and momentum con- 
servation laws but its evolution is inconsistent with gravity. This illustrates that in 
setting up the effective softening one should also ensure the consistency of com- 
puted forces with gravity. 

The uncertainty in the appropriate softening method may be resolved by resorting 
to a well-known physical interpretation of force softening. The appropriate phys- 
ical model is that simulation particles represent not point masses, but rather spa- 
tially extended clouds of finite density. This picture naturally leads to finite max- 
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imal pairwise forces : if p — p^ax towards the center of a cloud, then the force 
exerted on a test particle will reach a maximum and eventually vanish as the test 
particle approaches the interior and then the center of the cloud. Of course, to be 
self-consistent within this physical picture, we should not compute pairwise forces 
between particles treating one particle as an extended cloud and the other as a point 
particle, but rather treating both particles as extended. The pairwise interaction po- 
tential then becomes 



and the force is obtained by differentiation of the potential with respect to separa- 
tion between the centroids of clouds 1 and 2. 

Equation Q, while clearly the correct choice of force softening, may appear daunt- 
ing to evaluate. Pairwise force calculation is often the limiting step in N-body sim- 
ulations, and so an expression for the pairwise force that requires few floating point 
operations in its evaluation is an absolute requirement. The 6-dimensional inte- 
gral in Eqn. ^ would hardly appear promising in this regard, and so a simpler 
choice of softening, although less physically self-consistent, may seem more ap- 
pealing. However, we show in this paper that (perhaps surprisingly), the potential 
in Eqn. (O is in fact analytic for commonly used softening kernels. For the popular 
"spline" density kernel 0, the potential in Eqn. Q may be expressed with rational 
functions, while for Plummer force kernel the integral can be done numerically. 
This allows the efficient calculation of self-consistent forces between particles of 
unequal mass using our proposed W-shape and the extended Plummer force laws. 
Our force profiles can also be effectively used for particle splitting methods IIH ap- 
plied to pure gravity, or simulation of dark matter including pointlike objects such 
as super-massive black holes or stars. 

Section [2] gives general expressions for pairwise gravitational potential and force 
given the specified density shape for particles in the pair. A numeric integration 
is necessary for Plummer softening discussed in section [3l Closed form analytic 
solution for standard Wn shapes for n = 1,2,3,4 is presented in section |4l As 
n — i> oo, this solution asymptotically approaches gaussian softening discussed in 
section [5] In section [6] we compare all the discussed softening methods and relate 
them to each other. We conclude in section [8l 
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2 Review of Pairwise Forces 



Consider a simulation particle at position xq, of mass m and density profile p{x) = 
m W{x — xo), where the cloud shape W has unit normalization: 

W{x)d^x = l. (3) 



We will consider only spherically symmetric shapes, W{r) = W{r = \x — Xo\) 
to avoid generation of spin angular momenta from tidal torques. Then the Fourier 
transform of the shape similarly depends only on the magnitude k of the wavevec- 
tor k. 



W{k) = W{k) = J W{r)e-'''-^d^r 



(4) 



and is a real function. 

Given the particle's density profile, its gravitational potential is implicitly defined 
by 

VV(^) = 4vrp(f), (5) 



where we set Newton's gravitational constant Gat = 1 for simplicity. We can write 
an explicit expression for the potential using the Greens function for the Laplacian 
operator: 

(p{x) = / G^{x,y)p{y) d^y , (6) 



where the Greens function satisfies 

V'G^{x,y)=A7r6^'\x-y) . (7) 

Then, consistently with Eqn. we can write down the pairwise interaction energy 
between particles i and j as 



U.,{r-,) = U.,{r.,) = Jp.ix)^,{x)d'x_^ 

/-, h' 
W:{k)W,{k)G^{k)e^''^^j^^ 



2mimj 



TT 



W:{k)W,{k)Ukr,,)dk, 
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where we have used G^{k) = —Air/k'^, joix) = sm{x)/x, and rij = \fij\, fij = 

Differentiating with respect to the separation vector and flipping the sign gives 
pairwise force 

Fij(k) = -ikW;ik)Wjik)G^ik) , (9) 
applied on particle i by particle j. 

As an example, we can consider point particles, with W{f) = 5^^\r), or equiva- 
lentlyH^(A;) = 1. Then Eqn. ([8]) gives f/j^ = — m^mj/r or force Fjj = —mimjfij/rfj. 
More interesting applications will arise when we consider some other commonly 
used cloud profiles, discussed in the following sections. 



3 Plummer Softening 

Plummer softening Eqn. ([T]) leads to fourier component Fp(A;) = mimjATiike pKi{ke p) / k, 
where K^{x) is the modified Bessel function (using Eqns. 3.771.2 and 3.771.5 of 
M)- Comparing this with Eqn. Q we find W*{k)Wj{k) = kepKi{kep). 

This implies that Plummer force can be viewed as force between a cloud shape 
whose fourier component is Wi{k) = kepKi{kep) and a point mass Wj{k) = 1. 
Alternatively, Plummer force can be viewed as force between two identical Plum- 
mer cloud shapes 

Wp{ep, k) = ^kepKiikep) . (10) 

Figure [Hillustrates these alternatives in physical space perspective. 

Plummer force law can be consistently extended into pairs of particles on unequal 
softening by using Wp{ei, k) and Wp{e2-, k) in Eqn. ^ for computing potentials 
and forces for all pairs of particles. The integral leads to the Plummer force Eqn.© 
for pair of identical particles Si = 62 = sp. Numerical evaluation is required how- 
ever for pairs of unequal softenings. Our implementation in Appendix |7] includes 
potentials and forces for these Plummer shapes. 

As an alternative method of extending Plummer law into pairs of unequal softening, 
one would consider plugging a symmetric combination of Si and 82 into Eqn.dU). 
However, as we show in sections [5] and [6] such an approach leads to inconsistency 
of the resulting forces with gravity for pairs of unequal softenings. 
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Fig. 1. Force softening such as Plummer can be viewed as force between a cloud of shape 
W{k) > and a point mass (shown in physical space on top plate), or alternatively as force 
between two identical cloud shapes y'W{k) (bottom plate). We plot Plummer shapes with 
^yW{k) = Wp{ep, k) in this example, for ep = 1. 



4 Softening with Wn Clouds 



The assignment of particle mass to a regular grid is an important step in particle- 
mesh codes lHOl . and is achieved using one of several possible shapes: Nearest Grid 
Point (NGP), Cloud in Cell (CIC), Triangular-Shaped Cloud (TSC). In this section 
we discuss how these shapes are used to define our proposed W-shape softening 
and spline force softening previously used in literature. 



4.1 Hockney-Eastwood Cloud Shapes 



We can write the Hockney-Eastwood shapes in one dimension as 
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Scheme 


n 


Wnis),S > 
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1 
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CIC 


2 
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^ 8 ^3/2 (S) 


PQS 
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(|-^' + t)-^^i(«) 






+ 6 -^2(5) 



Table 1 

Hockney and Eastwood ifTOl cloud shapes in one dimension. Here, Hy{x) = Hq{x — y), 
where Hq is the Heaviside function defined as Ho{x) = for x < and Ho{x) = 1 for 
X > 0. 

The Wn{s = nx/h) clouds are characterized by two quantities: the scale length h, 
and the index n which controls the smoothness of the function. The function Wn{s) 
has n — \ continuous derivatives and disappears at x > 6/2. 

These Hockney-Eastwood cloud shapes are defined in one dimension (see Table[I]), 
however we can generalize them by replacing their argument s, a linear coordinate, 
with spherical radius r. Let us define 

Wr.{h,r) = —w^{nr/h) (12) 



where the prefactor is inserted to ensure that Wn is properly normalized. Because 
these cloud shapes have compact support, the force law they generate on a point 
test particle is exactly Newtonian for r > 6/2. 

Note that n = 4 corresponds to the so-called 'spline' softening used, for example, in 
smoothed particle hydrodynamics [|71 and pure gravity. Using b = 2has an exercise 
identically yields the density kernel in Eqn.(4) of [[TTl : using h = Ah yields Eqn.(ll) 
of [|T2l. 

To compute the interaction potential, we require an expression for the Fourier trans- 
form of Wn{b, r) to insert into Eqn. Using Eqns. ©, and we find 

Wn{h,k) = Wn{khl2n) (13) 



where 



w„M , llW f!ifi£V-' , (14) 



X \ X 
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and ji{x) = (sinx — xcosx)/x^ is the spherical Bessel function. As an exercise, 
using n = 2 and b = a one identically recovers Eqn.(A.16) of |fT3l . 



We have only verified this expression Eqn. (1141) for 1 < n < 4, however it may 
serve as our definition of the cloud shaped for an arbitrary n. Given this expres- 
sion for the smoothing kernel, the interaction potential for particles with smoothing 
scales bi and 62 becomes 

Un{biM,r)- - 



mirrij 

00 



IT J \ 2n ) \ 2n 





These integrals may be evaluated in closed form; the resulting expressions are 
lengthy and given in the appendix. For finite n, the force profile is Newtonian 
(/ oc r"^) outside r > (bi + b2)/2, and vanishes linearly (/ oc r) at the limit 
of zero separation. 

Interaction potential has a simpler form in the case of a pair consisting of two 
identical particles of smoothing scales b 

u^^\b,r) = unib,b,r) (16) 
or a pair including one point mass 




Note that as n — 00 at fixed b, Wnib, r) S^^\r). Therefore, even at fixed scale 
length b, the force profile approaches that of a point mass given this choice of 
normalization. 



4.2 W-shape Softening 



As noted above, the scale length b is not quite the softening length: the effective 
softening also depends upon the smoothing index n when we define the cloud pro- 
file to vanish exactly at r > 6/2. For convenience, we therefore choose to redefine 



^ Alternatively, one could define Wn(x) = ^ '^^ J which is consistent with taking a 
3-dimensional convolution. 
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the softening length to absorb this n dependence. We may do so by rescaling the 
softening length, writing 

b = Kf;iew (18) 



where the values of the coefficients K^^^^ for typical values n = \ — A are given in 
the table in Appendix IA.2[ The constant is chosen to make the interaction potential 
between two clouds of equal smoothing scale Sw at zero separation depend only 
on By/ and not on n, i.e. 

M(f)(6,r = 0) = — ^. (19) 



Figure [2] illustrates the density, potential, and force profiles for various n at fixed 
Sw, for pairs of identical particles. 

As may be apparent from the figure, as we take the limit n — >^ oo at fixed sw, the 
density profile converges to a gaussian discussed in section[5l Indeed, as suggested 
by table in Appendix I A . 2 1 the coefficient in Eqn. (fTSi) scales roughly as K^^ oc ^/n. 
Assuming this scaling for n — oo, using Eqns.dTSi). (fT3l) . (fT4l) . and taking the limit 
n — i> oo we arrive to the gaussian clouds Eqn. (I22l) . 

Next we consider the pairwise force between Wn shaped particles of unequal soft- 
ening length. Figure [3] shows examples of the force laws between particles of 
smoothing scales 8^ = 1 and Sw = l/o'; the different panels show different 
smoothing indices n. As the ratio g — > oo, this approaches the interaction between 
a Wn cloud and a point particle. As is apparent from the figure, the force profiles 
quickly converge to the asymptotic {n = oo) behaviors for both n > 1 and g > 1. 



4.3 Spline Softening 



As noted above, our n = 4 density shape corresponds to the spline kernels used 
in, e.g., SPH [I7I11I12II smoothing. However the spline force law does not exactly 
match ours. This is because the spline force used commonly in the literature does 
not correspond to the force law between two W4 clouds, but rather to the force law 
between a W4 cloud and a point particle. 

In spline softening, interparticle potential and force laws follow from Appendix |A.1.2[ 
Zero separation interaction potential becomes 

ui'\b,r = 0) = -^ = -f = -^, (20) 

6 50 ^spline 
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Fig. 2. Density, force and potential profiles for ew 
(ep = 1) profiles are shown. 

where we have used equation (IA.5I) and scaling 



1 and n = 1,2,3,4. Plummer 



'J -''■p4 ^spline • 



(21) 



It may first appear given the discussion in section [3] that spline softening can be 
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Fig. 3. Interparticle force laws in pair of particles of unequal softening lengths ew = 1 and 
sw = l/'Z- Solid lines correspond to different values of q = 1, . . . , 6 from bottom to top. 

viewed as force between the spherically symmetric clouds of shapes ^/W4. We can 
show however that spline softening is inconsistent with such model. Indeed, equa- 
tions (fT3l) and (fT4l) show W4 is negative for some k, hence ^^W^ is imaginary. On 
the other hand, equation dH) shows that the fourier component of any spherically 
symmetric cloud is real. By contradiction this proves that there is no possible so- 
lution for the density shape p(r) that for a pair of identical particles of this shape 
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would under gravitational interaction give us the force law used in the spline force 
softening. 

Whether or not this conclusion brings important consequences for simulations that 
use spline softening is debatable and is open for future tests. Using the extended 
Plummer, W-shape or gaussian softenings however allows one to immediately re- 
solve this uncertainty. 



5 Gaussian Softening 

Gaussian smoothing 



>Vg(^g, A;) =exp [-[kecf /'^t^] ■ (23) 

allows the simplest expression for softening between clouds of different smoothing 
scale and, as shown in section 14.21 is the n ^ oo limit of softening with our Wn 
shapes. For the interaction potential between two gaussian clouds of softenings ei 
and 82 whose centers are separated by distance r this gives 

Uciei, 62, r) = erf , (24) 

1" \ ^^sym / 



where 

^sym = ^sym(^l; ^2) = "i/ 77 (^1 + ^1) • (25) 




We have found that interaction potential in pair of gaussian clouds of softening 
scales £1 and £2 equals interaction potential between identical gaussian clouds each 
having softening scale esym- 

This simple prescription in Eqn. (l25l) however does not generalize to most other 
commonly used shapes. For a general cloud shape W{k) and general symmetric 
combination esyml^^i, £^2) the prescription is valid only if 

W^(£sym(^i,^2),fc) = [Wie,,k)W{e2,k)Y^' (26) 



for all k, as can be seen from Eqn. ([8]). The condition is exact for gaussian clouds 
and Esym given by Eqn. ( |25] ). but is not satisfied for other commonly used shapes: 
Plummer, spline or W-shapes. 
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6 Relations between Softenings 

In this section, we turn to the relation between our proposed W-shape (we assume 
n = 4), extended Plummer laws and the more familiar spline, Plummer and gaus- 
sian force softenings. 

How different, in practice, are our proposed profiles from previously used laws? 
To answer this, we first must normalize the various profiles to match each other as 
closely as possible. We do so by matching zero separation interaction potentials in 
pairs of identical particles. Normalized in this way 

£G = £P = = ^spline = £ , (27) 

which allows us to drop the subscripts. We see from figure |4] that our W^-shape 
profiles are quite close to spline and gaussian softening for the same e but diverge 
from Plummer profile. 

The close coincidence of W-shape and spline curves is consistent with the idea that 
our W-shapes reach gaussian in the limit n ^ oo. It may therefore appear first that 
the prescription e = \J {si + el)/^ of Eqn. (l25l) is the consistent way to generalize 
the interparticle laws to pairs of unequal softenings. As a measure of consistency, 
in figure[5]we plot the force laws found using this prescription against the "correct" 
force laws found in the result of double integration over particle shapes. The coin- 
cidence is identical for gaussian softening (not shown in figure), as we know from 
section [51 

For Plummer force law the solid lines are found by numerical integration using 
Eqns. ([8]) and (flOl) and the dashed lines are found by plugging Sp = ^ {ef + el)/^ 
into Eqn. ([T]). From the plot we find this prescription leads to up to 52% relative 
systematic inconsistency force. 

For W-shape softening this prescription results in at most 10% systematic increase 
over the self-consistent Poisson gravity force law. The latter is however easily com- 
putable using the expressions we provided, hence there is no advantage in using the 
simplification in the first place. 



7 Numerical Implementation 

The documented WSHAPE package (this paper uses version 1.0) is available at 
[http : / / www .graces . org/wshape, We provide numerical C-implementation 
of the potential and force laws for W-shape and extended Plummer. Also provided 
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Fig. 4. Plummer, W-shape, spline and gaussian force laws for e = 1. 
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Fig. 5. Force law in pairs of softenings ei = 1 and si = l/q,for q = 1, 2, oo in bottom 
to top order on each plate. Solid lines: found by double integration over particle shapes; 

Dashed lines: found by using the prescription e = \J (ef + 62) /2. The solid and dashed 
curves visibly overlap for q = 1. 
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within the package is the procedure used to arrive to expressions in Appendix |B] 
and values of coefficients in Table I A. ll 



8 Conclusions 

An N-body simulation is a simulation of massive particles under the influence of 
gravity. In cosmological simulation interparticle force is softened to avoid hard 
collisions. Simulation particles have often been viewed as particles of fixed density 
shapes in previous literature (e.g. expression for p(r) in Eqn. (5) of IfTTTl ). To ensure 
that the nature of interaction between particles in simulations remains gravitational 
we adpot this interpretation in which force softening corresponds to interaction 
between cloudlike particles of a fixed density shape. 

We prove that the Plummer interparticle force softening is consistent with this force 
softening model and generalize the Plummer force law into the case of unequal 
softening scales ei and 52- 

Interestingly, we mathematically prove that the previosuly used in literature (e.g. 
ifTTI ) method of spline softening is inconsistent with gravitational interaction in the 
sense that there is no possible solution for the density shape p{r) that for a pair of 
identical particles of this shape would under gravitational interaction give us the 
force law used for interparticle force law in spline softening. 

We provide the closed form solution for potential and force laws between widely 
known iy„-cloud shapes of general softenings ei and 62- 

Our generalized interparticle force laws should be useful for N-body simulations 
with adaptive mass refinement, in which particles of different mass interact gravita- 
tionally. Examples include simulations of the first collapsed objects in the universe 
m, the "Via Lactea" simulation [21, pure gravity extension of particle splitting (HI, 
and simulations of dark matter that includes pointlike objects, e.g. stars or super- 
massive black holes. 

As an easy way to extend the force law f{e, r) into pairs of particles of unequal 
softenings ei and £2 one would plug their symmetric combination of Si and 62 as 
the softening scale e. We found that using the simple gaussian motivated recipe e = 



y (^i + e2)/2 for simulations with adaptive mass resolution generally leads to sig- 
nificant systematic uncertainties for the pairs of particles of unequal softenings. The 
overall effect of these errors can be established by numerical convergence tests sim- 
ilar to [[T4|. However the problem is immediately resolved by using our proposed 
profiles. 

For convenience, we also made a numeric implementation of our analytic expres- 
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sions. The most efficient way to use tliese laws in an N-body simulation would be 
to pre-compute the interaction law once before running the simulation, and then to 
interpolate by look-up from the precomputed table. The GRACOS package IfTSl will 
implement these profiles for adaptive softening. 
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A Force and Potential Law in a Pair of Wn Cloud Shapes 

This section provides closed form expressions (for n = 1,2, 3, 4) for interaction 
potential laws in a pair of two Wn cloud shapes of scales bi and 62, a pair of two 
identical lV„-shape particles of scale b, and a pair consisting of an VF„-shape parti- 
cles and a point-like particle. Interparticle force laws follow immediately by differ- 
entiation with respect to separation r and flipping the sign. 

A.l Analytic Expression 

A. 1.1 Two Parametric Form 

Potential law Un{bi,b2,r) for interaction in a pair of two iy„-shape particles of 
scales bi and 62 and a fixed n is given by Eqn. (fTSi) . Its closed form solution is 
given as a finite sum over all integers i and j, with 

Un{bi,b2,r) = -^^A„,,(r) {r/a,r+^-\r/a2r^'-' (A.l) 

where as = bg/ (2M„) for s = 1, 2, and the values of positive integers M„ and An 
are given for each n in Appendix IB 

Coefficients Anij{r) depend on r through the Heaviside function Ho{r) via sum- 
mation over all integers p and q 

^nij{r) =Y.^ [Cnpq{i,j) + C„pg(j, i)] H{p, q, r) , (A.2) 

pg ^ 
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where H{p,q,r) = Hq^v — {pai + ^02)) for integers p and q, and Cnpqii,]) are 
constant non-zero integer coefficients, whose complete set for each n is presented in 
AppendixlHl The total number of terms in these expressions is finite. The potential 
is symmetric with respect to hi and 62, as should be. 



A. 1.2 One Parametric Forms 

Potential laws for interaction in a pair of two identical iy„-shape particles of scale h 
(case c = 2), or one PF„-shape particle of scale h and a point-like particle (case 
c = 1) are found in Eqn. (fTTl ) and (fT6l ). Their closed form solutions are given by 

«W(6,r) = --^E^2(r) (r/a) ("+-1)^+2- , (A.3) 
An r i 



where a = 6/(2 Af„). The values of positive integers M„ and A^^^ are given for 
each n and c = 1, 2 in AppendixlHl Coefficients (r) depend on r through the 
Heaviside function Ho{r) via 

A^hr) = EcS(^H{p,r) (A.4) 

p 



where H{p, r) = Ho{r — pa) for an integer p, and C^^{i) are constant non-zero 
integer coefficients, whose complete set for each n is also presented in AppendixlHl 

As an exercise compute uf\4h, r), using the above expressions and table values. 
You should arrive to the right hand side of Eqn.(A2) of [|T2ll . which is the potential 
kernel in the traditional spline softening described in Section |431 



A.2 Asymptotic Expressions 



Potential laws in section lATTI are pure Newtonian beyond separation distance when 
particle last overlap. Indeed 625 = — l/ratr > (6i+62)/2, and r) = 

— l/ratr>c6/2. On the other hand, for small separations we have analytically 
fore = 1,2 



nW(6,r) = -i^W/6 
/W(6,r) = -rfr/63 



as r <C & , 



(A.5) 



where numerical values of coefficients are given in Tables lA.ll 
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Table A.l 

Table of coefficients for Eqn. (IA.5l ) 



Pair: one W-shape and one point particle (c = 1) 



n 


1 


2 


3 


4 




3 

8 


4 
32 


39/8 
54 


28/5 
256/3 



Pair of identical W-shape particles (c = 2) 



n 


1 


2 


3 


4 


J^^pn 
^(2) 


12/5 
8 


104/35 

64/5 


124/35 
774/35 


70016/17325 
31424/945 



Figure lA.ll shows the density, potential and force laws for n = 1 — 4 for pair of 
identical particles of smoothing scale 6=1. Note the different amplitudes for the 
potential and linear force law at small r for different n, which is consistent with 
growing values of K!^^ and Kf^ for n = 1, 2, 3, 4. 



B Tables of Coefficients 



This subsection provides tables of coefficients for analytic expressions for poten- 
tials in Section IATI for n = 1, 2, 3, 4. The procedure used to find these table values 
is given in Appendix |7l 



B. 1 Spherical Top-hat - Shaped Particles (Wi -Shape ) 

Values n = 1, M„ = 1 are applicable for the entire subsection lB.li 



B.1.1 Point-like Test Particle 

4^) = 2 
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n = 


4 






















/ ^~ 


3 \ 










' / n = 


2 \ 










n = 


1 


V 
















-0.4 


-0.2 C 


) 0.2 


0.4 0. 




n=4 

-.!)=2 
-.n=1. 










Plummer, e=i 



















-0.5 

logio('-) 



0.5 



logio('-) 



Fig. A.l. Density, force and potential laws for 6 = 1 for unsealed W^-shapes for 
n = 1, 2, 3, 4. Plummer law for e = 1 is also shown. 

The following table lists values of coefficients C'>^^{%)\ the numbers, labeling the 
rows and columns denote p and %. 



cm 



np 

2 3 




1 



-1 3 
1 -3 2 
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B.1.2 Identical Particles 

= 160 



The following table lists values of coefficients C^^^{i)\ the numbers, labeling the 
rows and columns denote p and i. 



cm 








2 3 


5 


6 





-1 


30 -80 


192 




2 


1 


-30 80 


-192 


160 



B.1.3 General Case 
An = 160 

The following two tables list values of coefficients Cnoo{i,j) and Cnu{i,j)', the 
numbers, labeling the rows and columns denote i and j. 



CnOoihj) 








2 3 4 5 


6 





1 


-30 -80 -90 -48 


-10 


2 




90 240 90 




3 




-80 

















2 3 4 5 


6 





1 


-30 80 -90 48 


-10 


2 




90 -240 90 




3 




80 





The rest of the non zero coefficients in Eqn. (IA.2I) are given by the following rela- 
tions 
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B.2 Cone - Shaped Particles (W2-Shape) 



Values n = 2, Mn = 1 are applicable for the entire subsection |B.2[ 



B.2.1 Point-like Test Particle 

A^:^ = 1 



The following table lists values of coefficients C^J{i); the numbers labeling the 
rows and columns denote p and i. 





np 


(0 









1 


3 4 





1 


-2 


2 


1 


-1 


2 


-2 1 



B.2. 2 Identical Particles 

^(,2) = 140 



The following table lists values of coefficients C^J{i); the numbers, labeling the 
rows and columns in the following table, denote p and i respectively. 



Cnp (?) 








1 


2 


3 


4 


5 


6 


7 


8 





3 


-8 


-14 


56 




-112 




208 




1 


-4 


16 




-112 


280 


-336 


224 


-80 


12 


2 


1 


-8 


14 


56 


-280 


448 


-224 


-128 


128 



B.2. 3 General Case 
An = 140 
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The following three tables list values of coefficients Cnooihj)^ Cnm{i,j), andC„ii(i, j); 
the numbers, labeling the rows and columns denote i and j respectively. 



CnOoihj) 





1 


3 4 5 


6 7 8 





1 -8 


56 140 168 


112 40 6 


1 


14 


— 140 —280 —252 


— 112 —20 


3 




140 280 84 




4 




-70 








CnOl{i,j) 






1 


3 4 5 


6 7 8 





-2 8 


-56 140 -168 


112 -40 6 






Cnu{i,j) 






1 


3 4 5 


6 7 8 





1 -8 


56 -140 168 - 


-112 40 -6 


1 


14 - 


-140 280 -252 


112 -20 


3 




140 -280 84 




4 




70 





The rest of the non zero coefficients are given by the following relations. 

Cnlo(j,0 = C'nOl(^,j) 
C„,i,_i(z,j) = {-ly Cnll{^,J) 

B.3 TSC- Shape Particles (W^-Shape) 

Values n = 3, M„ = 3 are applicable for the entire subsection lB.3[ 

B.3.1 Point-like Test Particle 
= 160 
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The following table lists values of coefficients C^J{i); the numbers, labeling the 
rows and columns denote p and i respectively. 

12 4 5 

2 -20 130 

1 -3 10 -10 5 -2 
3 1 -10 30 -135 162 

B.3.2 Identical Particles 
^(f) = 12902400 

Table lB.ll lists values of coefficients C^^{i). 

B.3.3 General Case 
An = 12902400 

Table |B?2] lists values of coefficients C«ii(z, j). The rest of the non zero coeffi- 
cients Cnpq{i, j) are given by the following relations 



Cn,i,-i{hj) = i-iy^' Cnu{^,J), 



C„33(^,j) = M-2+*+^- all(^,J), 



B.4 Cubic Spline - Shaped Particles (W4 -Shape) 

Values n = 4, M„ = 2 are applicable for the entire subsection IB .41 
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B.4. 1 Point-like Test Particle 
= 30 

The following table lists values of coefficients C^l^{i)\ the numbers, labeling the 
rows and columns denote p and i. 








1 


2 


3 


5 


6 





-3 


9 




-20 


42 




1 


4 


-18 


30 


-20 


6 


-2 


2 


-1 


9 


-30 


40 


-48 


32 



B.4.2 Identical Particles 
^^2) = 1663200 

Table lB.3l lists values of coefficients C^J{i). 

B.4.3 General Case 
An = 1663200 



Table lists values of coefficients Cnoo{i,j), Cnoi{i,j) and Cnu{i,j)- The rest 
of the coefficients Cnpq{i, j) are given by the following relations: 



Cnlo(j,0 = 




CnOl 


%3) 


Cn, -1,1 




= i-iy 


Cnll 


%J) 


Cn02(«,j) = 


_J\^-2+i+j 

n 


CnOl 


%J) 


Cn,2,-2 


%J) 


= M-^+'+^{- 


Cnll 


[hj) 


Cn20(j,^) = 




CnOl 


%J) 


Cn, -2,2 


%J) 


= M-^+*+^'(- 


Cnll 


[hj) 


Cnuihj) = 


n 


Cnll 


%J) 


Cn, 1,-2 


%J) 




-ly Cnll 




Cn2lii,j) = 


n 


Cnll 


%J) 


Cn, -2,1 


%J) 




-ly Cnll 




Cn22ii,j) = 


n 


Cnll 


%J) 


Cn,2,-1 


%J) 




ly Cnll 




Cn,l-l(«,j) = 




Cnll 


%J) 


Cn, -1,2 


%J) 




-ly Cnll 
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Table B.l 

Table of coefficients Cnp{i), where n = 3. Numbers labeling rows and columns in the following 
table denote p and i. 








1 


2 


3 


4 


5 


6 


7 


8 


9 10 





-10 


80 


180 


-2880 




48384 




-660480 




7618560 


2 


15 


-200 


810 


1440 


-26880 


120960 


-302400 


468480 


-449280 


245760 -58880 


4 


-6 


160 


-1620 


5760 


26880 


-387072 


1935360 


-5406720 


8847360 


-7864320 2883584 


6 


1 


-40 


630 


-4320 




217728 


-1632960 


5598720 


-8398080 


10077696 



Table B.2 

Table of coefficients Cnoo{hj) and C„ii(i,j), where n = 3. Numbers, labeling the rows and 
columns in the following two tables, denote i and j. 








2 


4 


5 


6 


7 


8 9 10 





4 


-360 


21840 


161280 


609840 


1397760 


1967400 1574400 551096 


2 




5040 


-327600 


-1612800 


-3659040 


-4193280 


-1967400 


4 






2129400 


10483200 


7927920 






5 








-6451200 









Cnll{i,j) 








1 


2 


4 


5 


6 


7 


8 


9 


10 





9 


-120 


270 


-1260 


3024 


-3780 


2880 


-1350 


360 


-42 


1 




360 


-1440 


5040 


-10080 


10080 


-5760 


1800 


-240 




2 






1260 


-6300 


10080 


-7560 


2880 


-450 






4 








3150 


-5040 


1260 










5 










1008 
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Table B.3 

Table of coefficients Cnp{i), where n = 4. The numbers, labeling the rows and columns in the 
following table, denote p and i. 








1 


2 


3 4 


5 


6 7 


8 9 


10 


11 12 





35 


-180 


-165 


2200 


-15840 


95040 


-432080 




1680384 


1 


-56 


504 


-1716 


1760 5940 


-26928 


53592 -66528 


55440 -31240 


11484 


-2496 244 


2 


28 


-504 


3828 - 


-14960 23760 


50688 - 


-384384 1064448 


-1774080 1914880 


-1317888 


528384 -94208 


3 


-8 


216 


-2508 15840 -53460 42768 


449064 -2309472 5773680 -8660520 7794468 


-3779136 708588 


4 


1 


-36 


561 


-4840 23760 


-50688- 


-118272 1216512 


-4055040 7208960 


-6488064 1572864 1048576 



Table B.4 

Table of coefficients Cnoo(^> j), CnOiihj) and C„ii(i,i), where n = 4. The numbers, labeling 
the rows and columns in the following three tables, denote i and j. 





1 


3 


5 


6 


7 


8 


9 


10 


11 


12 





9 -108 


1320 


-33264 


-166320 


-441936 


-748440 


-838200 


-605880 


-257544 


-49104 


1 


297 


-5940 


116424 


498960 


1104840 


1496880 


1257300 


605880 


128772 




3 




18480 


-388080 


-1108800 


-1473120 


-997920 


-279400 








5 






814968 


2328480 


1031184 












6 








-831600 















CnOlihj) 








1 


2 


3 


5 


6 


7 8 


9 10 


11 


12 





-24 


216 


-792 


1320 


-4752 


11088 


-14256 11880 


-6600 2376 


-504 


48 



Cnll{i,j) 








1 


2 


3 


5 


6 


7 


8 


9 


10 


11 


12 





16 


-288 


1056 


-1760 


6336 


-14784 


19008 


-15840 


8800 


-3168 


672 


-64 


1 




1188 


-7920 


11880 


-33264 


66528 


-71280 


47520 


-19800 


4752 


-504 




2 






11880 


-31680 


66528 


-110880 


95040 


-47520 


13200 


-1584 






3 








18480 


-55440 


73920 


-47520 


15840 


-2200 








5 










16632 


-22176 


4752 












6 












3696 
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